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Abstract 

Nearly fifty years after the introduction of skew Brownian motion by Ito and 
McKean (1963), the first passage time distribution remains unknown. In this 
paper, we generalize results of Pitman and Yor (2001) and Csaki and Hu (2004) 
to derive formulae for the distribution of ranked excursion heights of skew 
Brownian motion. We then derive the first passage time distribution as a 
simple corollary. 
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1. Introduction 

In this paper, we obtain formulae for the distributions of first passage time and ranked excursion 
heights of skew Brownian motion. Since Ito and McKean (1963) first introduced skew Brownian 
motion, numerous papers in the mathematics literature have highlighted the special structure of the 
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process and connected it to diverse applications; see, for example, Walsh (1978), Lc Gall (1984), 
Barlow ct al. (1989), Oukninc (1990), Decamps et al. (2006), Ramirez et al. (2006), Ramirez et al. 
(2008), and Ramirez (2010). 

Skew Brownian motion is broadly applicable to diffusion problems in which the diffusion coeffi- 
cient is discontinuous in space. For example, Appuhamillage et al. (2009) analyzed mathematical 
properties of skew Brownian motion as they relate to the convection-dispersion movement of solute 
through porous media in the presence of a sharp interface. Their work was motivated by experiments 
in hydrology demonstrating a fundamental asymmetry in the breakthrough curves of solute crossing 
a sharp interface in opposite directions (Berkowitz et al., 2009). In this context, the first passage 
time distribution of skew Brownian motion describes the breakthrough times of solute injected on 
one side of an interface and observed on the other side (Appuhamillage et al., 2010). 

Similar situations arise problems in physical and natural sciences. For example, Brownian motion 
has been widely used, and critiqued, as a model of animal movement (Turchin, 1996), (Blackwcll, 
1997). However, there is an emerging consensus among ecologists that landscape heterogeneity is 
a necessary ingredient of movement models (Dalziel et al., 2008). In particular, sharp interfaces 
delimiting two different movement regimes occur at the boundaries of habitat patches (Schultz and 
Crone, 2001; Turchin, 1991; Wiens and Milne, 1989), or at other environmental discontinuities such 
as those in ocean temperatures (Pinaud ct al., 2005) or the level of surface chlorophyll in oceans 
(Polovina et al., 2001). In an ecological context, McKcnzie ct al. (2009) advocate the use of first 
passage time to model the time required for an organism to first reach a specified location in the 
landscape (see also Fauchald and Tveraa (2003)), and give as an example the foraging behavior of 
a predator searching for stationary prey items. 

The main contribution of the present paper is to derive the first passage time distribution of 
skew Brownian motion. We achieve this by first characterizing the distribution of ranked excursion 
heights of skew Brownian motion. This result generalizes formuale of Pitman and Yor (2001) for the 
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distribution of ranked excursion heights of the standard Brownian bridge, and analogous formulae 
presented by Csaki and Hu (2004) for Brownian motion. We then apply our formulae for the ranked 
excursion heights of skew Brownian motion to derive our main result on first passage time. 

The paper is organized as follows. In Section 2, we state the main results. In Section 3, we 
develop a coupled construction for two different skew Brownian motion processes with different skew 
parameters that leads to an important relationship between distributions of ranked excursion heights 
of the two processes, stated in Theorem 3.1. In Section 4, we prove the main results as corollaries 
to Theorem 3.1. In Section 5, we present several examples that demonstrate the calculations that 
are possible using the first passage time density, and also the assymetry inherent in the first passage 
time of particles crossing an interface in opposite directions. 



2. Preliminaries and Main Results 

To set some notation and basic definitions, let B — {B t : t > 0} be the standard Brownian 
motion process on a probability space (f2, T, P) and let Ji, J2, . . . denote the excursion intervals of 
the reflected process {\B t \ : t > 0}. For a G (0, 1), let {A^m ' m = 0, 1, . . .} be a sequence of i.i.d. 
±1 Bernoulli random variables with P(A$ = 1) = a. Define the a-skew Brownian motion process 
B^ started at by 

OO 

B[ a) = J2 l Jm (t)A^\B t \, (1) 

m— 1 

where Is denotes the indicator function of the set S. 



Now let 

M[ a \t) > M[ a \t) > ... > 

be the ranked decreasing sequence of excursion heights sup se j- mn r 0jt ] Bg ranging over all m such 
that J m n [0,t] is nonempty. Note that a negative excursion has height zero, and that the height of 
the final excursion is included in the ranked list even if that excursion is incomplete. Our first main 
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result gives the distribution of ranked excursion heights. 



Theorem 2.1. Fix y > and t > 0. Then for each j — 1, 2, the distribution of Af- (i) is given 



by the formula 



Po(M ( ] a \t)>y) = 
where <&(•) is the standard normal distribution function. 
Now let 

= inf{s > : flW = j/} 

denote the first time for a-skew Brownian motion to reach y and let (x, ?/, t) denote the first 
passage time density to y at time t of a-skew Brownian motion started at x. When a = 1/2, this 
is the well known first passage time density f(x,y,t) for Brownian motion (e.g., see page 30 of 
Bhattacharya and Waymire (2009)): 



) = £ 2 ( ■ _ !) (i - 2a ) h ~ 3 (i - *((2* - i)y/Vt)) , 

h=l ^ ' 



\y - gj 

2^t 3 /2 



exp {- ~ 2^} 



Notice that f ( - a \x,y,t)dt = P x (T^ a) € di). Our second main result gives formulae for the first 
passage time density. 

Theorem 2.2. Fix t > 0. TTien 



/or x < < y 



f {a) {x,y,t) = i 



9™ (*) + W ~ £ — siu{ ^ X)n } g W * /or 0<x<y 

71=1 ^ 

f(x,y,t) for < y < x 



and 
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/or y < < x 



n=i 71-71 y 

f(-x, -y,t) for x < y < 



h («) o VVi o V7-1 k ~ ( 2 J ~ f (x - (2j - l)y) 2 \ 
where g x ,y — 2a y (1 — 2aj J — _ — — — exp| — — j /or x < t/ie quantity 

K>n,y(t) is the density at time t of the exponential distribution with parameter X(n,y) — ™ 2 ™ 2 ; and 

. , 7r f n 2 n 2 t] . ( 7r(y — x)n] 

h XtV (t) — — nexp< — — > sim > is the well known formula (e.g. see page 296 in 

V n=1 2y y 

Feller (1968)) for the probability that Brownian motion started at x reaches zero before reaching y, 
and that this event occurs in the time interval dt. 



The cases < y < x and x < y < in the first passage time density are clear because all paths 
starting at x reach y without hitting zero and hence they are all Brownian motion paths. In cases, 
x < < y and y < < x, all paths must cross zero and densities are obtained as a straightforward 
corollary of Theorem 2.1 and the fact that P a {T^ a) e dt) = P (T { ^y a) e dt). The situation is most 
complicated when < x < y and y < x < 0, and one must consider two types of paths from x to y: 
those that cross zero, and those that reach y before they reach zero. 

Notice that when a — 1/2, we recover existing results for standard Brownian motion. Namely, 
from Theorem 2.1, we recover the distribution of ranked excursion heights stated in Theorem 3.1 of 
Csaki and Hu (2004), and from Theorem 2.2 we recover the well known first passage time distribution 
of standard Brownian motion (this fact is not immediately obvious, but nonetheless true, in most 
complicated cases when < x < y and y < x < 0). 
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3. Relating excursion heights for and 



Let < a < j3 < 1. Consider the following coupled construction of a-skew and /3-skew Brownian 
motion. Let B be the standard Brownian motion process and let — {Am^ : m = 0, 1, . . .} be 
independently chosen excursion signs so that Equation (1) yields an instance of /3-skew Brownian 
motion. 

Next, let {A^ /f5) : m = 0, 1, . . .} be a sequence of i.i.d. ±1 Bernoulli random variables independent 
of A<#> and B with P(Aln /P) = 1) = a//3. Define A!$ as follows 



4(a) 



1 A iP) - 1 A (a/P) - 1 

r "ra — r,^im — -L- 



-1 otherwise. 



By construction, the sequence {if : m = 0, 1, . . .} consists of i.i.d. ±1 Bernoulli random variables 
that are independent of B with P{Am ^ = 1) = a. Hence, by using the variables i™' as the excursion 
signs in Equation (1), we obtain an instance B^ of a-skew Brownian motion. 

We think of this as a two-step process: first, construct B^ by independently setting each excur- 
sion of \B\ to be positive with probability (3; then, for each positive excursion of B^\ independently 
decide whether to keep it positive (with probability a//3), or flip it to be negative (with probability 
1 - a/13). 

The following theorem is motivated by this coupled construction. 

Theorem 3.1. Fix y > 0, t > and a, /3 € (0,1). For each j = 1,2,..., the following relation 
between ranked excursion heights of a- and /3-skew Brownian motions holds. 

P (M^(t) > y) = £ ( h : \ ) (l - f ) ~ J (f) J Po(Mf (t) > y). (2) 

Before giving the proof of Theorem 3.1, we state the following lemma from Pitman and Yor 
(2001), as we use it in the proof. 
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Lemma 3.1. (Pitman and Yor (2001), Lemma 9.) Let 



h = ^2 f , \a m , k = 0, 1, 

m=0 ^ ' 



be the binomial moments of a nonnegative sequence (a m , m = 0, 1, ...). Let B(9) := X^fclo ^fe^ and 
suppose B{6{) < oo for some Q\ > 1. Then 

am = E(-l) fc - m Q)bfe, m = 0,l,..., 
where the series is absolutely convergent. 

Proof of Theorem 3.1. For a < ft, we have by the two-step construction of the excursion sign 
im' that Mj(t) — Mg}(t), where Hj has a negative binomial distribution: 



(;•:;)(- frw- 



Hence 



Po(M^(i) > y) = f] J) (l - I)" J '(f ) J P (M(« (t) > y). (3) 



For ft < a, the relation can be inverted by an application of Lemma 3.1. Let k := j — 1 and 
m := h — 1. Then one can write (2) as 

Po(Ml%(t) >y)=jt (*?) (l - 1 ) m ~\f ) k+1 Po(M^ +1 (t) > y). (4) 

m=0 ^ ' 

We then apply Lemma 3.1 to the sequences 

h := (l - f) fc (f)" fc_1 P (M£ ) 1 (t) > y), a m := (l - ^ P^M^t) > y). 
After simplifying, we obtain 

p o( Mf ( t) > s) = £ (j: 1 ) (i - ^^(fy^'w > „). (5) 
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4. Proofs of Main Theorems 

We now observe that the main results announced in the introduction will follow as corollaries to 
Theorem 3.1. We first prove the Theorem 2.1. 

Proof of Theorem 2.1. By Theorem 3.1 in Csaki and Hu (2004), 

P Q {Mf' 2 \t) > y) = 2(l - *((2j - 1)^)) (6) 

The result is immediate from Theorem 3.1 by taking j3 — 1/2 in (2). 

We now use Theorem 2.1 and the following corollary to compute the distribution of first passage 
time asserted in Theorem 2.2. 

Corollary 4.1. Fix t > 0. Then 

f2aE^ 1 (l-2^- 1 ^cx P {- ^_nv)l} dt fory>0 

P (T^edt)=l J (7) 

(2(1 - a) £r=i(2a - exp{- ^^^} dt fory<0 

Proof. For the case y > and t > 0, we have the following relation between the distributions of 
Ty 01 ^ and the highest excursion of skew Brownian motion started at 0: 

Pq(tW <t) = P (M[ a \t)>y). 

Thus using Theorem 2.1, one has 

P {T^ <t) = P Q {M[ a \t)>y) 

= 4aV(l-2a) /l - 1 / -= expl \ dz. 

rt=l Vt 

The result is immediate after taking the derivative of the above expression with respect to t. 
For the case y < 0, use the relation Po(Ty € dt) = FoC^ly G dt) and the case y > 0. 

Proof of Theorem 2.2. Let T a = Ty 1 ^ denote the first time for standard Brownian motion to 

\y\ ( y 1 i 

reach y. Recalling that Po{T y € dt) = — ^=-^-^expj — — j dt; e.g. see page 30 in Bhattacharya 
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and Waymirc (2009), one can write equation (7) as 

2a££Li(l " 2a) /l - 1 P (7 1 (2/l _ 1)y G dt) for y > 

P (T^edt)=l (8) 

[2(1 - a) Er=i(2« ~ l) fe - 1 ^o(T(2/ l -i) 3/ G dt) for y < 
Now note that Tg is distributed as To under P x for x ^ 0, < a < 1. So clearly for t > 0, one 



has 



P E (T (a) > t) = P x (To > t). 



(9) 



(10) 



We prove the case y > in the Theorem 2.2. The case y < is similar. 
Case x < < y : 

Using the strong Markov property of skew Brownian motion, 

P X (T^ > t) = / P^To > t - s)P (tW g ds). 
Jo 

Then from the first case of equation (8), one has 

00 -x 

P x {T< a ) > t) = 2«V(1 - 2a)' 1 - 1 / P X {T Q > t - s)P (T {2h _ 1)y G ds) 
h=i Ja 

OO 

= 2a ^(1 - 2a)' l - 1 P,(T (2 , l _ 1)!/ > t) 

h=l 

By differentiating the above expression with respect to i and recalling 
P X (T V G dt) = I^T X [ exp(- ( y ~ x } X dt, one has 

G dt) = 2.£(1 - 2^- l l ^ v ^; /2 1)yl ex P {- ( - (2/ ;-^ )2 } dt (11) 

1 

Case < x < y : 
Observe that 

P,(TW G dt) = P X (T^ G dt, (T (Q) < t)) + P,(T^) G dt, (T (Q) > t)) (12) 

We state the following formula for P x (Tq G dt, (T y > t)); e.g. see page 296 in Feller (1968), as we 
use it to compute (12) 

P K (r edt,(T y >t)) = ^2jncxp[- — ^-jsinj j dt (13) 

y n= i 2/ 2/ 
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Since the skew Brownian motion is Brownian motion until it reaches zero for the first time and 
from the reflection principle of Brownian motion, one can write the second term of the right hand 
side of equation (12) using equation (13) as follows: 

P x {T { y a) G dt, (T (Q) >t)) = P x (T y G dt, (T > t)) 

= P y - x {T edt,{T y >t)) (14) 

= 7L nffl p{- -^} sm l — — I dt 

For the first term of the right hand side of equation (12), notice that 



P x (Tt ] < < t) 



E, 



1 lT^ ) <Ti c " ) <t] 



E, 



E, 



-.(a) 



L [T o w <TW<t]l x i - L [T (Q, <T^ ) ] 



[Tr'<Tr'<*] 



t (q) - s 1 



+ / E a 

Jo 



1 [T (Q) <T^ ) <t]l T ) - S ' 1 [T ( < Q) >T^ ) ] 



P,(T (Q) ed s ,(r (Q) >tW)) 



Using the strong Markov property of skew Brownian motion and since E x 



[T (Q, <T ! $ a, <i]l J -0 



-,(«) _ 



= 0, one has 



P X {T^ ] < tM < f) 



1 [T <°) <t _ s] 



Tx(T (Q) Gd S ,(T (Q) <TW)) 



P (T( Q ) < t - s)P x (Tt ] G ds, (T (Q) < TW)) 



Using equation (8) and again from the fact that the skew Brownian motion is Brownian motion 
until it reaches zero for the first time, and P x (Tq <T") = P X (T < T y ) (note here < x < y), one 
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P X (T M < <t)= 2aV(l - 2a) h - 1 P Q (T {2h _ 1)y < t - s)P x (T € ds, (T < T y )) 
Jo h=i 

oo „t 

= 2aV(l-2af- 1 / P (T (2/l _ 1)s/ < t - s)P x (T E ds, (T < T y )) 

u-1 JO 



h=l 

oo 



oo „ t 

2a^(l-2a)' 1 - 1 / P {T (2h _ 1)y < t - s) P X {T e ds) 

h 1 JO 

(15) 



OO 



- 2a V(l - 2a)' 1 - 1 / P (T (2h _ 1)3 , < « - s)P x (T e ds, (T > T„)) 

OO 

2a - 2a)' l - 1 P ;E (T (2/l _ 1)y < i) 

oo „ t 

- 2a V(l - 2a)' 1 - 1 / P {T (2h _ 1)y < t - s)P x (T Q G ds, (T > T„)) 
fc=l • y ° 



For /i > 1, the convolution integral in the sum of the second term in the above equation can be 
written using (13) as; 



•t 

Po(T (2 h-i)y < t - s)P x (T Q € ds, (T > T y )) 

i 

= / P x (To<t-s,(T >T v ))P (T i2h _ 1)y eds) 
Jo 

Jo nn 1 V J 

' 7T71 L W J v 

n=l y 

E2 r 7r 2 72 2 t s \ ( Tr(y — x)n l /"* r7r 2 n 2 si ,„ , , 

^} Sm { J V^} X exp{^}Po(T (2fe _ 1)y e rfs) 



Then one has 
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P X {T^ < < t) 

OO 

= 2a ^(1 - 2a) h - 1 P x (T (2h _ 1)y < t) 

h=l 

oo oo { \ 

00 00 n 2 2-1 / \ 

h=ln=l y y 

(n 2 n 2 s^\ 

x X exp {"2^r o(T ( 2h - i )i' eds) 



By differentiating the above equation with respect to t, one has 



2a ^(1 - 2a)' l - 1 P :r (r (2 ^ 1) , e <ft) 



oo oo 



2a££(l-2a)'^sin{ 

h=ln=l 



2 r 7r(w — x)n 
sm<| 

7rn 



}p°(7< 



7r 2 n 2 £i . r7r(?/ — a;)n 



oo oo 



2y 2 



2/1-1)9 G d*) 

} 



(18) 



x exp< 



2/y 



h=\ n=l 



- 2a££(l-2a)'--exp{-_ }sm{^-^} 



c* , r 2„2. 

x / exp-j 



Recalling that P x (T y e d£) = X [ exp(- ^-—^ } (it, and by (12), (14), (18), one has 

V27T t 3 " L 2t ) 
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P X {T^ e dt) 

Al — 1 

OO OO 9 9, / \ 

„ ~ xh -i^n ( 7T n t ~] . ( niy — xm} 

- 2aEE( 1 - 2 «) h -V eXP {-^} Sm { J V^} 

ft,=ln=l a y f 

/•* f7r 2 n 2 s-| (2ft- l)y ( ((2h-\)y) 2 } , , 

Jo I 2y 2 J ^s 3 / 2 I 2s J 



OO 



7T v-^ f 7T n t] . f 7T( 7/ — acini . 

7^{ nexp \ ^} sm l y 1 dt 



2aV(l-2a)^^^=iMe X p{- (» " ^ " W \ 
^ \/2^t 3 / 2 I 2i J 



, , /2¥t 3 / 2 

Al — 1 



/•*7r 2 n 2 r 7r 2 7i 2 (t-sV» „ Nh i(2/i-l)i/ f 2/i-l y 2 ! , , 

x / ^exp^ K — '- } 2a > (1 - 2a) h ~ 1 ^^ — ^ exp<^ - ^ \ ds dt 

J 2y 2 y \ 2y 2 I ^ ' y/2ir s 3 /2 v \ 2 t I 



7r ^— v r 7r n ti . f7r(y — a;)ni , 

+ ^ nexp r ^l sm l — „ — I rfi 

(19) 

Case < y < x : 

Notice that, in this case all skew Brownian motion paths till the first passage time to y are away 
from 0. Thus we have 

\y~ x \ „_r (y-x) 2 , dt 



P (T& G dt) - lV ~ Xl exnl- [V ~ X) \ 



5. Asymmetries in First Passage Time Density 

A basic property of the first passage time density of regular Brownian motion is that it is 
symmetric in x and y, i.e., that P x (T y G dt) = P y (T x G dt) for any x,y G K and t > 0. A 
practically important feature of skew Brownian motion is that it introduces an asymmetry in the 
first passage time density for x and y on opposite sides of the origin. Appuhamillagc et al. (2010) 
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first demonstrated this by proving that there is a stochastic ordering between the random variables 
with densities P- y (Ty a ^ £ dt) and P y {T_ y ' £ dt) when a ^ 1/2. In simple terms, it takes longer 
to cross from negative to positive than to cross from positive to negative when a < 1/2, and the 
opposite is true when a > 1/2. 

Our results quantify these relationships further by explicitly giving the densities in each case. We 
conclude by illustrating the numerically computed densities P_i(T-l € dt) and PiiT^ e dt) for 
different values of a. In Figure 5(a) we see that P-i(T^ a) E dt) < Px(T^ G dt) for a < 1/2. In 
Figure 5(b), we recover the symmetry of Brownian motion when a = 1/2. In Figure 5(c), we see 
that P-i(Tl a) e dt) > Pi(T^ e dt) for a > 1/2. 

a = 0.1 a = 0.25 




(c) a > 1/2 

FIGURE 1: The densities P-i(T^ a) G dt) (solid) and Pi(T^ G dt) (dashed) for different values of a. 
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